One-dimensional Kondo lattice model as a Tomonaga-Luttinger liquid 
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Arguments are presented that in the one-dimensional Kondo lattice model /-electron spins partici- 
pate in filling of the Fermi sea. It is shown that in its paramagnetic phase this model belongs to the 
spin-1/2 Tomonaga-Luttinger liquid universality class. The ratio of the spin and charge velocities 
Va/vp and K p are estimated to be of the order of (Tk/cf) 1 ■ 
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I. INTRODUCTION 

One can use a one dimensional Kondo Lattice (KL) 
model as a toy model to study the long standing prob- 
lem of whether localized electrons determine the volume 
of the Fermi surface. The numerical results pf, [0 show 
that in the region of the phase diagram where KL model 
belongs to the universality class of spin-1/2 Tomonaga- 
Luttinger (TL) liquid [p[ the Friedel oscillations are char- 
acterized by the large Fermi vector. So it seems that the 
/-electrons do participate in the Fermi surface formation. 

In this paper we undertake a further study of the TL 
phase of the KL model. Let the reader recall that the 
spin-1/2 TL liquid critical point is characterized by two 
parameters: the ratio of spin to charge density wave ve- 
locities v a jv p and the number K p which parameterizes 
scaling dimensions in the charge sector (the similar pa- 
rameter in the spin sector is fixed by the SU(2) symmetry, 
K a = 1). The dynamical spin and charge susceptibilities 
at low [uj, q) are given by 



X*(u n ,q) = - 



7T q 2 v a + wl/vo- 



x P (un,q) 



2K n 



ir q 2 v p + ujI/v p 



(1) 



(2) 



Thus if we shall manage to find these two parameters, the 
characterization of the low energy sector of KL model is 
complete. 



II. LUTTINGER LIQUID PARAMETERS 
OBTAINED BY THE 1/JV EXPANSION 

The Hamiltonian of the one-dimensional KL model is 
F=-t^(aJ CT a i+1(T + H.c.) + J^-Sf^ (3) 



where a ifJ (a^) is the creation (annihilation) opera- 
tor of a conduction electron at the ?th site, and sf = 

( l / 2 )T,aa> a L T aa' a «?'' with T aa[ (/* = X i V' Z ) bein S the 

Pauli matrices, are the spin density operators of the con- 
duction electrons. The spin densities are coupled to the 
localized spins S% through an antiferromagnetic exchange 
coupling J. 

In order to obtain analytical results we shall extend the 
symmetry of the KL model to the SU(N) and resort to 
the 1 /TV-expansion (see Q and ||). The corresponding 
Lagrangian density is 



C = a*[d T + e(x)]a j + f*d T f j + 

iA(/;/;-«iV)-^(a*/,-)(/ fc *a fc ) 



(4) 



Here the dynamical field A(r, x) is introduced to en- 
force the local constraint of the fermion occupation num- 
ber. The number q remains finite when N — * oo. 

Next we decouple the interaction by the Hubbard- 
Stratonovich transformation: 



J^( a *jfj)(fk a k) 



N rz +V{a * fj)+V * {f * aj) . 



(5) 



The resulting partition function is gauge invariant: 

V(T,x)^V(T,x)e- [,KT ' x) , 
A(t, x) — v A(t, x) — 8 t (J)(t, x) . 



(6) 



It is convenient to choose the gauge where the field V is 
real. We choose the following parametrizations: 



V(t,x) = V yJl + [r(T,x)/V VN\ , 
iA = T K + iu/VN . 



(7) 
(8) 



where Vb is the saddle point value of V which we shall 
determine later and r is a new field chosen in such a way 
that its measure of integration is trivial. 

We shall expand the partition function around its sad- 
dle point: 



V = Vb , iA = T K . 

Expanding to the second order in r we get 



(9) 



£■ — -Co + C-int 

C a = a*[d T + e{x)] aj + f*(d r + T K )f 3 + V^f, + c.c.) 

(10) 



C 



int ~ AJ' r 2VN VUjJJ 



T 17/ 

(: </ J+ c.c.:) + -=:/;/, 



(11) 



where the dots mean that the average is substructed: 
: A := A — (A) and the r 2 term comes from the expan- 
sion of the square root in the expression for V (pi) . The 
saddle point parameters Vb and Tk are determined self- 
consistently by vanishing of the terms linear in r and u: 



N 



^2J<tf(»)/i(n)> = «, 



1 N 

-Y 



H(n)fj(n)) + (//(n)oi(n))] 



2Vb 
J 



(12) 



(13) 



In the leading order in 1/JV the spectrum is determined 
by the saddle point. This gives us a great advantage 
because the saddle point describes the large Fermi sur- 
face. The single electron spectrum has the following well- 
familiar form: 



E±(p) = [e(p) + T K ]/2 ± ^/[e(p)-T K ] 2 /4 + V 2 ■ (14) 

Substituting the saddle point Green's functions into Eqs. 
(H H) we get 

g = P (0)V o 2 /T K , T K = £>exp[-l/p(0)J] (15) 

where p(0) is the bare density of states per one channel, 
D is the bandwidth. The expression for the new Fermi 
vector is 

k F = k ( p ] + irq . (16) 

It follows from this equation that the charge susceptibil- 
ity remains unaffected by the presence of the spins. 

Only the mode E-(p) crosses the chemical potential. 
Near the Fermi points the spectrum can be linearized: 

E_(p)^±v*(pTPF), v* = v FP (0)T K q . (17) 

Now one can calculate the spin- and charge density re- 
sponse functions (g) directly. The calculation of the spin- 
spin correlation function is straightforward: in the lead- 
ing order in 1/N the only contribution comes from the 



polarization loop of two Gf = ((// + ))-functions. The 
result reproduces Eq. (fy) with v a = v*. 

Calculation of the charge response function is more 
complicated. We chose the following approach: first we 
shall integrate out the high energy degrees of freedom in 
the partition function and obtain the effective action for 
the low energy sector; then we shall bosonize this action 
and obtain the parameters of the TL liquid. To do the in- 
tegration it is convenient to diagonilize the saddle point 
Hamiltonian and to express the fermionic operators in 
terms of the new annihilation operators A±Ak) corre- 
sponding to excitations with the dispersion E±(k) (since 
the transformation is diagonal in the flavour indices, we 
shall omit them): 



a(k) = ^c k A + {k) + y/0 k A_(k) 
f(k) = -^p k A + (k) + V^ k A-(k) 



(18) 



where 



a k + f3 k = 1, p k = - < 1- 



2 I j[e{k)-T K \*+W* 



(19) 



Substituting these expressions into Eq. ( |ll| ) and omitting 
the terms containing only A + we get: 



dxCi 



Li + L 2 



(20) 



*i=E 



r(-q)r(q) 
13 



k,q 



2VN 



2m{q)yJ{3 k+q a k [A% Ak + q)A_ tj (k) + c.c] (21) 



L? 






[\/a F l3 F r(q) + ia F u(q)} 



(22) 



where a F , p F are taken at the Fermi surface: a F ~ 
1, I3 F ~ T K /D. 

Integrating over A+ we get in the leading order in 1/N 
the following action for the fields r and u: 



S ^ff = «E n ^ 9 ) 



[4r(-w, -q)r{uj, q) + u(-w, -q)u(u>, q)] . (23) 

To get the effective action for the low-lying excitations 
we need to know the function T1(uj, q) for the area around 
q = and for q = 2k F . The result is 11(0, 0) = p(0). 

We bosonize the fermionic operators: 



^2A*_ !J (k + q)A_ tj (k) = iy/N/^q$ p (q) (\q\ « k F ) 

(24) 

where $ p is the charge field and integrate over u and 
r. Since f3 F is so small, the largest contribution to the 
effective action comes from the fluctuations of the w-field. 
The bosonized version of the effective action in the charge 
sector is given by 



S eff = jdTdxl^(d T $t 



1 



271-/9(0) 



(0**p) 2 ]. (25) 



From this action one can derive the canonical expres- 
sion for the charge susceptibility and K p . At least in the 
leading order in 1/N the result does not depend on N: 



v p = y/v*/irp(0) 



and 



K p = ^/np(0)v* 



(26) 



(27) 



III. DENSITY MATRIX RENORMALIZATION 
GROUP STUDY 

In order to check the validity of the large N results for 
the N = 2 case we numerically estimate the TL liquid pa- 
rameters making use of the density matrix renormaliza- 
tion group (DMRG) || . This method is the most suitable 
for studying long range and low-energy properties since 
it allows one to study long chains iteratively enlarging 
system size and to obtain the ground state wave func- 
tion with only small systematic errors, which can be es- 
timated from the eigenvalues of the density matrix. The 
obtained results are consistent with the above arguments 
and indicate K p <C 1/2 in the Kondo limit. 

Now we shall describe results of the numerical anal- 
ysis of the model (||). The paramagnetic metallic state 
of this model, which is expected to be a TL liquid, is 
realized only in the region of rather weak exchange cou- 
pling away from both the half-filling (n c = 1) and the 
low carrier density limit (n c — > 0). The ground state is 
always insulating at half-filling and ferromagnetic both 
in the strong coupling limit (J — > oo) for general carrier 
densities (n c 7^ 1) and in the low carrier density limit M. 

We first calculate spin excitation gap A s and difference 
of chemical potentials p+ — pi- as a function of the system 
size L. As expected, both A s and p + — /i_ ( Fig. 1 (a) 
and (b) for the case of n c — 2/3 and J = 1.8i, 2.0t) 
vanish in the bulk limit (L — > 00), which confirms that 
the paramagnetic phase of the KLM is a TL liquid. 

The finite size corrections of /i + — fi- and A s in 
Fig. 1 are related to the charge susceptibility and the 
spin velocity, respectively. Since we have used open 
boundary conditions A a (L) — v a Ak(L) = v a Tr/L, and 



fi+(L) - ii-(L) = An c (L)/ X p = 2/( Xp L). The ob- 
taincd values are shown in Table I. Once we have ob- 
tained v G then we can calculate \<j through the relation 
K a = 7rw (T x CT /2, see Eq. (|l|). Because the SU(2) symme- 
try in the spin space guarantees K a = 1 , rather large \„ 
is obtained as is shown in Table I. This large x<j is nat- 
urally expected because there are macroscopic number 
of almost free spins in both weak and strong coupling 
regions. The / spins are almost but not exactly inde- 
pendent with each other: in the weak coupling region, L 
almost free / spins, and in the strong coupling region, 
i(l — n c ) /-spins unpaired with conduction electrons. 

Now we discuss the charge susceptibility. In the strong 
coupling limit it tends to the value for the free spinless 
fermions; x" 1 = 7ri sin (7T — 7rn c ) . On the other hand 
in the weak coupling limit we expect a J independent 
charge susceptibility as is predicted by the Gutzwiller 
type variational calculations S; Xp 1 — Trfsin (7rn c /2). 
The density n c = 2/3 is rather special in the sense that 
the values expected for the strong coupling limit and the 
weak coupling limit are the same. Thus we expect % p de- 
pends only weakly on J. In general, in the weak coupling 
limit we have an asymptotic form of the charge velocity 
which is proportional to K p as 



= 2K p t sin (7rn c /2) 



(28) 



from the relation K p — nVpXp/2, Eq. (pi). However, we 
have to be careful close to the half-filling where the charge 
susceptibility tends to diverge owing to the charge gap at 
half-filling. 

The estimation of the correlation exponent is one of the 
most difficult calculations even by the DMRG method. 
In order to estimate K p we need to see long range be- 
haviors of the system with sufficient accuracy. In the 
present study we use asymptotic form of the Friedel os- 
cillations because they are numerically more reliable than 
long range off-diagonal correlations. 

The Friedel oscillations are density oscillations induced 
by a local perturbation. In a TL liquid, power low 
anomalies in correlation functions naturally reflect them- 
selves in the Friedel oscillations; the Friedel oscillations 
induced by an impurity potential are 

Sp(x) - d cos{2k F x)x^ l ~ K ^/ 2 + C 2 cos{Ak F x)x- 2K » 

(29) 

as a function of the distance x from the impurity M O] . 
and analogously, spin density oscillations induced by a 
local magnetic field behave as 



<j(x) ~ D\ cos(2k F x)x 



-k 



(30) 



Thus, we can determine K p from the asymptotic form of 
the oscillations. 

Fig. 2 shows induced charge- and spin-density Friedel 
oscillations of the KLM obtained by the DMRG for 
J = 2.5i at n c = 6/7. The Fourier components of 
spin-density Friedel oscillations for J = 1.8i, 2.5i, 2.5i at 



n c = 2/3,4/5,6/7, respectively, are also shown in Fig. 3. 
The charge density Friedel oscillations are induced natu- 
rally by the open boundary conditions of the system and 
the spin density oscillations are introduced by applying 
local magnetic fields at the both ends. As is already 
shown for n c = 4/5 in the previous work |l[], the pe- 
riod of the oscillations are explained by the assumption 
of the spin-1/2 TL liquid with the large Fermi surface, 
kp = 7r(l +n c )/2, which includes / spin densities as well 
as the density of conduction electrons. 

Now we calculate the correlation exponent, K p . In 
order to obtain K p , we simply use the slope of the enve- 
lope function of the charge density oscillations assuming 
that its decay is proportional to x~ 2Kp , because domi- 
nant component of the oscillations is the 4&f component 
even for the case of J = 1.54. In Fig. 4, the obtained 
K p for the exchange coupling from J = 4.04 to 1.54 at 
n c = 2/3 are presented. Since the 2kp spin density oscil- 
lations decay much slower than the charge density oscil- 
lations, it is not possible to determine K p from the spin 
density oscillations in the present system size. However, 
the slower decay of the spin density oscillations is consis- 
tent with the TL liquid prediction, Eq. (pCf), which gives 
smaller exponent, x~ Kp . 

As is clearly seen in Fig. 4, K p is always smaller than 
1/2 and monotonically decreases with decreasing J. In 
the strong coupling limit, the conduction electrons and 
the localized / spins form local singlets leading to a com- 
plete spin-charge separation. Since the charge part is 
described by the free spinless fermions, K p = 1/2 is ob- 
tained in the strong coupling limit as in the case of the 
infinite-[7 Hubbard model. With decreasing J from in- 
finity the repulsive interaction between the neghboring 
spinless fermions is introduced in the leading order of 
t/J. Thus the situation is similar to the large-I/ Hub- 
bard model with nearest neighbor repulsions whose K p 
is smaller than 1/2 |UJ]|. 

In Fig. 4 we find a small discontinuity at J = 2.44. 
This is due to the phase transition from the ferromagnetic 
state to the paramagnetic one. Since this transition is of 
the first order accompanied by a jump in the total spin 
quantum number, from S = L(l — n c )/2 to 0, or 1/2 
with L being the number of the sites, it is natural that 
the K p also shows a jump at the critical value J c . In 
order to confirm the discontinuity we have calculated the 
K p in both ferromagnetic and the paramagnetic states 
at J = 2.44 which is near but smaller than the critical 
point. The K p in the ferromagnetic state is calculated by 
setting the total S z being L(l — n c )/2 which is the total 
spin in the ferromagnetic state. 

In contrast to the slow decrease of K p above the critical 
J c , a rather sharp decrease is observed below J c , and 
the Kp becomes smaller than 1/3 which means that the 
long range behavior of the charge-charge correlation is 
governed by the Akp oscillations. The dominance of the 
Akp oscillations is a characteristic feature of this new 
class of spin-1/2 TL liquid. 

With further decreasing J, the K p seems to cross the 



value 3 — 2\/2 ~ 0.17. Since the exponent of the power 
low anomaly in the momentum distribution function, a, 
is given by a — (K p +l/K p —2)/A, the power low anomaly 
is removed below this point and we cannot see clear Fermi 
surface any more. It is very difficult to observe clear 
Friedel oscillations for smaller K p than 0.17. 



IV. CONCLUSIONS 

In conclusion we have established that in the area of 
phase diagram where the one dimensional Kondo lattice 
is paramagnetic, it belongs to the universality class of 
spin-1/2 TL liquids. The /-electrons do take part in for- 
mation of the Fermi surface. According to the Luttinger 
theorem the volume of the Fermi sea is determined by 
those branches of the spectrum which cross the chemi- 
cal potential. Despite of the fact that most of the spec- 
tral weight of the /-electrons is concentrated far from 
the chemical potential, they do have access to it via the 
Kondo resonance. We bring attention of the reader to 
the fact that the Luttinger theorem does not require ex- 
istence of a pole in the single electron Green's function 
and therefore can be applied outside of the Fermi liquid 
domain. In particular, the system under consideration 
belongs to the spin-1/2 TL liquid universality class. It is 
a rather peculiar member of this class since K p is small. 
How small is not entirely clear; the analytical calculations 
give K p — Vfj/vp (see Eqs.rt2fl, |27|)) which is one order of 
magnitude smaller than the values obtained numerically 
(see Table I). This may be due to the inaccuracy of the 
1/iV-approximation; it is more likely however, that the 
maximal system size available for the numerical calcula- 
tions is not big enough to penetrate to the asymptotic 
region. Thus the numerical values of K p given in the Ta- 
ble should be considered as upper limits. This may ap- 
pear unusual to those who consider the Hubbard model 
as a typical example of TL liquid. The smallness of K p 
clearly originates from the nonlocality of the effective in- 
teractions in space and time. In this sense the KL model 
is similar to Charge Density Wave systems where the in- 
teractions are also retarded being carried by low energy 
optical phonons. In these systems K p <C 1 14 1, |lq| . 
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TABLE I. Luttinger liquid parameters of the one dimen- 
sional Kondo lattice model. The carrier density n c is 2/3. 
The energy unit is t. The errors are estimated from the am- 
biguity of the power low decay of the charge density Friedel 
oscillations. 
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FIG. 1. (a) Size dependence of the difference of the chemi- 
cal potentials, (j,+ — fi-, in the one dimensional Kondo lattice 
model. 2/u+(L) = E g (n c = n° + 2/L,L) - E g (n c = n°,L). 



2fi- (L) = E g (n c = n°,L)- E g (n c 



2/L,L). E g 



,L) 



is the ground state energy at the carrier density n c in the sys- 



= 2/3. (b) Size dependence of the spin 
= l,L)-E g {Sl ot = 0,L). E g (Sl ot ,L) 



tem of length L. n c 

gap A S (L) = E g (Sf 

is the lowest energy in the Hilbert space of total spin S*°*. 

n c = 2/3. The energy unit is t. Typical truncation errors in 

the DMRG calculations are 10 -4 . 



FIG. 2. (a) Charge density Friedel oscillations induced by 
the open boundary conditions. The system size is 70 sites, 
(b) Spin density Friedel oscillations induced by applying lo- 
cal magnetic fields at the both ends. The strength of the 
local magnetic fields is 0.2*. Typical truncation errors in the 
DMRG calculations are 1 x 10~ 6 for J = 2.5*. 



FIG. 3. Fourier components of the spin density Friedel 
oscillations. 



FIG. 4. Correlation exponent K p estimated from the de- 
cay rate of the charge density Friedel oscillations. The er- 
rorbars are determined from the ambiguity of the power low 
fitting. n c = 2/3. J is in units of *. 
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